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Abstract. We present the results of the asteroseismological analysis of two rich DAVs, G38-29 and R808, recent targets 
of the Whole Earth Telescope. 20 periods between 413 s and 1089 s were found in G38-29's pulsation spectrum, while 
R808 is an even richer pulsator, with 24 periods between 404 s and 1 144 s. Traditionally, DAVs that have been analyzed 
asteroseismologically have had less than half a dozen modes. Such a large number of modes presents a special challenge 
to white dwarf asteroseismology, but at the same time has the potential to yield a detailed picture of the interior chemical 
make-up of DAVs. We explore this possibility by varying the core profiles as well as the layer masses. We use an iterative grid 
search approach to find best fit models for G38-29 and R808 and comment on some of the intricacies of fine grid searches in 
white dwarf asteroseismology. 

Keywords: White Dwarfs, Asteroseismology, G38-29, R808 
PACS: 95.75.-z; 97.20.Rp 



ASTROPHYSICAL CONTEXT 

G38-29 and R808 are high amplitude, cool hydrogen 
dominated atmosphere pulsating white dwarfs (cDAVs). 
cDAVs are of interest for studying convection using the 
shape of their highly non-linear light curves [1]. In addi- 
tion, G38-29 and R808 are rich white dwarf pulsators, 
with over 20 periods present in their power spectrum. 
This makes them ideal subjects for asteroseismological 
studies. With the exception of hot pulsating white dwarfs 
[e.g. 2] , such rich pulsating white dwarfs have not been 
analyzed asteroseismologically before. Because of the 
large number of periods, asteroseismological identifica- 
tion of the modes in these stars is necessary in order to 
keep the study of convection using non-linear light curve 
fitting techniques computationally tractable. 

With such a large number of modes, we can begin to 
probe the interior structure of white dwarfs in more de- 
tail. Cool pulsating white dwarfs have the advantage of 
being simpler to model than hot white dwarfs. The chem- 
ical elements in their interiors are nearly fully settled and 
much of their interiors close to fully degenerate. In addi- 
tion, there are over 100 known DAVs [3] as opposed to 
5 known DOVs [4]. Not all DAVs are rich pulsators, but 
with improved observational techniques, there is a poten- 
tial to increase the number of modes observed in these 
stars, providing useful data for detailed asteroseismolog- 
ical analyses of white dwarf interiors. 

In this paper, we present preliminary asteroseismolog- 
ical analyses of G38-29 and R808, based on Whole Earth 
Telescope (WET) campaigns performed in fall 2007 and 
spring 2008 [5]. We are able to derive stellar parame- 
ters such as mass and effective temperature for each star 
as well as internal structure parameters. We suggest an 



identification for the modes in each star. 



SPECTROSCOPY AND POWER 
SPECTRA 

According to spectroscopy [6], G38-29 has a tempera- 
ture of 11,180 K and a logg of 7.91, which translates to 
a mass of 0.55 Mq. It was the object of a small WET 
campaign in November 2007. Analysis of the data re- 
vealed 20 independent modes, with periods ranging be- 
tween 413.307 s and 1089.39 s. Of special interest in the 
power spectrum is an ^ = 1 triplet, identified as such from 
average period spacing arguments. The triplet is split by 
7 jUHz, leading to a rotational period of 21 hrs for the 
star and providing ground for the identification of addi- 
tional multiplets (see Table 2). Out of the 20 independent 
modes, 15 are candidate m=0 modes. 

R808 is similar in temperature to G38-29 (11,160 K), 
but more massive with a logg of 8.04, corresponding to a 
mass of 0.63 Mq. It was the object of a WET campaign 
in April 2008, where 25 independent modes were found 
in its somewhat noisier period spectrum. We tentatively 
identified A t = 2 multiplets and I £ = I multiplet (see 
Table 2), consistent with a rotation period of 18 hours. 
The rotational sphtting analysis leaves 18 possible m=0 
modes. 



ASTEROSEISMOLOGICAL ANALYSIS 

The periods observed in white dwarfs are g-modes, 
where the restoring force is buoyancy. In a completely 
homogeneous star, g-modes would be evenly spaced in 



period. In reality, white dwarf's are differentiated and 
the chemical transition zones induce a departure of the 
modes from their even spacings. The amount by which 
the modes deviate from their even spacing provides clues 
to the interior chemical structure of white dwarfs. In the 
presence of rotation, the modes get evenly split in fre- 
quency. The magnitude of the frequency split depends 
both on the rotational period and on the i identification 
of the mode [7]. 

A rotational frequency splitting of 7 /jHz corresponds 
to a Am = 1 period split of 1 second for a 400 s mode and 
8 seconds for a 1100 s mode. For 1=2 modes, the corre- 
sponding frequency split is near 1 1 /iHz, corresponding 
to 2 and 13 seconds spUt in periods for the 400 s and 
1100 s modes respectively. These values are to compare 
with an average period spacing of ~42 s for (=\ modes 
and ~27 s for 1=2 modes. For asteroseismological fits, 
this means that lower radial overtone modes (lower k) are 
less sensitive to the exact m identification than higher k 
modes. Fitting lower k modes first therefore minimizes 
modeUng uncertainties due to our lack of knowledge of 
what the m identification of the modes is. On the other 
hand, lower k modes are also more sensitive to core struc- 
ture in the model and magnify our a priori ignorance of 
the core chemical profiles. As a first pass, we assumed 
that the observed modes were m=0 modes and looked 
for an asymptotic period spacing in the observed period 
spectra of G38-29 and R808. 

For our asteroseismological fits, we used the White 
Dwarf Evolution code (WDEC) to generate white dwarf 
models. The WDEC is described in detail in Lamb and 
van Horn [8] and Wood [9], and more recent modifica- 
tions in Bischoff-Kim et al. [10]. We used an iterative 
grid search method varying 6 parameters, listed in Ta- 
ble 1, starting with a low resolution grid covering a wide 
region of parameter space in a range of effective temper- 
ature and mass suggested by the spectroscopic values. 
The grid is presented in Table 1. Based on the fits found, 
we successively refined the grid, focusing on the regions 
of parameter space where the best fit models resided at 
each iteration. The final grids have a resolution of 50 K 
in effective temperature, 0.002 in mass for G38-29, 
and 0.005 for R808, 0.002 (in the log) for the heUum 
and hydrogen layer masses, 0.02 in central oxygen abun- 
dance and 0.02 for Xfm, the point where the oxygen 
abundance first starts to drop down. 



at mass point 0.68 M*. For R808, the corresponding pa- 
rameters are 11,250 K, 0.675 Mq, lO-^-^s, lO'^-^^^ 
0.10, andZ/„ = 0.84. 

The quality of each period fit for G38-29 is illustrated 
in Figure 1 . For (. = \ modes, the uncertain m identifica- 
tions can lead to uncertainties in model parameters, but 
not in I and k identification. For t = 2 modes, the situa- 
tion becomes more ambiguous for higher k modes, as the 
frequency splitting translates to a period splitting that is 
of the same order as the average period spacing. Consec- 
utive multiplets start to overlap. 
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FIGURE 1. Quality of fit for the G38-29. The "error bars" 
are not error bars, but indicate the extent of the Am = ±1 
rotational period splitting mentioned in the text. Note that this 
splitting is significantly less than the period spacing for 1= \ 
modes, though they can lead to significant uncertainties in the 
model periods and parameters if misidentified in m value. For 
the £ = 2 modes, the frequency splitting translates to a period 
splitting that is of the same order as the average period spacing. 
The situation is aggravated by the fact that m = ±2 modes are 
allowed as well. 



RESULTS 



DISCUSSION 



The results of the fits for G38-29 and R808 are presented 
in Table 2. For G38-29, the best fit model has a temper- 
ature of 11,550 K, a mass of 0.642 M©, a heUum layer 
mass of 10~^-^^, a hydrogen layer mass of lO""*-^^. The 
central oxygen abundance is 0.80 and it starts dropping 



In grid searches, the question often arises "how fine 
is fine enough?". How finely do we need to sample 
parameter space in order to find any minimum present? 
SampUng tests are needed to fully answer that question, 



TABLE 1. Parameters for the starting grid 



Description 


Parameter 


Range 


Step size 


Effective temperature 


Teff 


10,600 -11,800 K 


200 K 


Stellar mass 


M* 


0.500 - 0.700 Mq 


0.010 Mq 


Helium layer mass 


MHe 


10-2.00 . 10-2-40 


0.20 in the log 


Hydrogen layer mass 


Mh 


down to 10- ■'■00 


0.20 in the log 


Central oxygen abundance 


Xo 


0.00-1.00 


0.10 


Point where oxygen abundance 








first drops to zero 


Qftn 


0.10-0.80 


0.10 



TABLE 2. Mode identification and best fit model periods for G38-29 and 
R808 

G38-29 R808 



Oh *j f rvp H 


Mndel 


I 




ni 


Oh*^p^vpH 


Mndel 






HI 


Period 


Period 








Period 


Period 








432 


428 


1 


7 





405 


406 


1 


6 





545 




1 


10 


+1 


745 


741 


1 


14 





547 


545 


1 


10 





875 


870 


1 


17 





549 




1 


10 


-1 


912 


912 


1 


18 





706 


705 


1 


14 





915 




1 


18 


-1 


709 




1 


14 


-1 












962 


957 


1 


20 





952 


954 


1 


19 





1002 


1002 




21 





1040 


11)42 


1 


21 





1082 






23 















1089 


1086 




23 


+1 












413 


409 


2 


14 





511 


514 


2 


17 















629 




2 


22 


+1 












632 


637 


2 


22 















796 


788 


2 


28 















843 


837 


2 


30 















860 


862 


2 


31 





840 


844 


2 


32 





878 




2 


32 


+1 












899 




2 


32 


-2 












908 




2 


33 


+1 












916 


914 


2 


33 















923 




2 


33 


-1 


900 


900 


2 


34 





952 




2 


35 


+1 


923 


927 


2 


35 





961 


965 


2 


35 





945 


940 


2 


36 















964 


963 


2 


37 





1011 


1016 


2 


37 





980 




2 


38 


+1 












990 


992 


2 


38 





1042 


1042 


2 


38 





1016 


1022 


2 


39 





1067 


1067 


2 


39 















1091 


1092 


2 


40 





1087 


1087 


2 


42 





1144 


1144 


2 


42 






though we can check that our solution approaches the 
minimum monotonically. The step sizes for the final, 
high resolution grid quoted in the previous section appear 
sufficiently small. This also means that if one wanted 
to create a grid that covers the parameter ranges listed 
in Table 1, one would have to calculate ~ 10^*^ models. 
This is not manageable and the recourse is either to use 
"smarter" search algorithms such as genetic algorithms 
[11] or successively zooming in on the best fit. The 
two methods are complementary to each other. Genetic 
algorithms, given enough modes, are effective at finding 



a global minimum while grid searches help refining the 
minimum and provide more flexibility in trying different 
mode identifications. 

The results presented here are preliminary, as we 
have made the assumption that all observed modes that 
were not members of multiplets were m=0 modes. As 
noted earlier, this assumption does not significantly in- 
fluence the fits for k<10 modes, but becomes impor- 
tant for higher k modes. Observationally, m=0 modes 
are not necessarily the higher amplitude modes in mul- 
tiplets suggesting that single modes may not be m=0 



modes. Lightcurve fitting techniques provide an indepen- 
dent way of mode identification [12]. It would be worth 
trying more general fits and using light curve fitting tech- 
niques for G38-29, where the observed frequencies are 
better determined. 
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